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pi • Abstract. We present a study on the spectra for heavy flavor (charm and 

D \ bottom) decayed electrons in 200 GeV p+p collisions and provide the relative 

, ' contributions of charm and bottom hadrons from the PYTHIA calculations. 

The results suggest that the crossing point of the electron spectra from charm 
and bottom decays is above 7 GeV / c and bottom contribution is not dominant 
\ for electron < 3 GeV/c. The upper limit of the relative cross section ratio is 

reported as cTf,^ /dec < (0.49±0.09±0.09)%. We also compare the W2 distribution 
\ from simulation to the experimental data in 200 GeV Au+Au collisions and 

^SJ ' estimate the possible charm V2. 
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^ , 1. Introduction 

H , 

Cu ^ Due to the absence of the measurement of B-mesons and precise measurement of 

D-mesons, it is difficult to separate bottom and charm contributions experimentally 
in current non-photonic electron measurements for both spectra and elliptic flow V2 ■ 
As discussed previously, the suppression behavior of heavy quarks is quite different 
from light quarks due to the "dead cone" effect [[T], and this is especially true for 
the bottom quark. Even when the elastic energy loss is included, the bottom quark 
still loses much less energy. The bottom contribution may reduce the energy loss 
of non-photonic electrons from heavy flavor decays. But recent measurements show 
that the suppression of the non-photonic electron Raa is as large as light hadrons [ 
[12] . Both the theoretical result with charm energy loss only and the theoretical 
calculations with charm+bottom energy loss by assuming large q or counting elastic 
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energy loss can describe the data within errors [ [51 [31 HI [31 [B] . 

Recently, PHENIX has measured the non-photonic electron V2 [ [7] . The ob- 
served large elliptic flow of the non-photonic electron may indicate strong coupling 
of heavy quarks with medium. There are many theoretical calculations for the 
non- photonic electron V2 , such as charm thermal-f flow model [ [8j , A multi-phase 
transition (AMPT) model assume cross section cjp—lO mb [[9], resonance states 
of D-/B- mesons [[lO], etc. The comparison with theories also showes that both 
the model results with charm only and the results with charm+bottom have good 
agreement with data within errors. 

Thus, the puzzle of the bottom contributions in non-photonic electron spectra 
and V2 still remains. We present the following method to estimate the bottom 
contributions and to study the possible charm V2- 

2. Fit to non-photonic electron spectrum and relative cross 
section ratio 

The non-photonic electron spectrum up to 10 GeV/c has been measured by STAR 
experiment in 200 GeV p+p collisions. The idea is that we use the sum of electron 
spectra from both charm and bottom decays in PYTHIA model [ [IT] to fit the 
STAR p+p data [ [T^ to extract the fraction of the bottom contribution. Since 
the D-mesons and their decay electrons spectra from default PYTHIA parameters 
are soft [ [T3] , a modified Peterson Fragment Function (FF) and the high pT tuned 
parameter are used to make spectra harder to be comparable with the form factor 
decays [[T4]. 

Table [H lists the parameter initialization for PYTHIA 6.131: 



Table 1. PYTHIA parameters for heavy flavor decays. 



Parameter 


Value 


MSEL 


4 (charm), 5 (bottom) 


quark mass 


= 1.25, mb = 4.8 (GeV) 


parton dist. function 


CTEQ5L 


Q'^ scale 


4 


K factor 


3.5 


{Kt) 


1.5 


Peterson Frag, function 


e = IQ-"" 


high PT timed PARP(67) 
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Fig.[T](a) shows ihepx distributions of the heavy flavor hadrons and their decay 
electrons from PYTHIA with above parameters. The D-meson spectrum, shown as 
the hatched band, is normalized to 



dN/dy = dN/dy{D°)/{NMn)/RdAu/R, 



(1) 
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Fig. 1. Panel (a): D-/B- mesons and their decay electron spectra from PYTHIA. 
The B + 13 e fit to STAR non-photonic electron data in p+p collisions. Panel 
(b): The relative spectra ratio, upper limit of _B — > e contributions as a function of 
Pt- 



where dN/dy{D°) = 0.028 ± 0.004 ± 0.008 measured in rf+Au collisions [ ITS]. 
(Nbin) = 7.5 ± 0.4 in d+Au collisions. RdAu = 1-3 ± 0.3 [ 16 . R factor stands 
for fraction in total charmed hadrons, the fragmentation ratio R{c D^) = 
Njjo/Ncc = 0.54 ± 0.05 [ 17 . All these normalization errors are propagated into 
the uncertainty band of the D-meson spectrum. The curve in this band is the lower 
limit of the D-meson spectrum in our simulation. Correspondingly, its decay elec- 
tron spectrum is shown as the solid band. The non-photonic electron spectrum 
measured in p+p collisions at STAR [ I12j is shown as the open squares. The decay 
electron band alone can describe the data, indicating that the contribution of elec- 
trons from bottom decay could be very small. In order to estimate the upper limit of 
bottom contribution, we use the lower limit of the decay electron spectrum, shown 
as the open circles. B-meson spectrum (solid curve) and its decay electron spec- 
trum (open triangles) are normalized by varying the ratio of cjjji/iJcc- The summed 
spectrum (solid circles) by combining the lower limit oi D ^ e and i? ^ e is used 
to fit STAR data in p+p collisions, and then the upper limit of _B ^ e contribution 
will be extracted. 

Fig. [2] (a) shows the fit as a function of the unique variable (j^ijocc- The best 
fit with a minimum jndf = 16.6/14 gives the upper limit of the total cross section 
ratio as u^ijucc < (0.49 ± 0.09 ± 0.09)%. The first term of the errors is calculated 
from — Xmin + 1- The second term is from the 15% normalization error of 
the dN/dy converted to total cross sections due to the uncertainties of the model 
dependent rapidity distributions [ 16J. Fig. [2](b) shows the B-/D- mesons rapidity 
distributions from PYTHIA. The cross section ratio from FONLL calculation is 
0.18%-2.6% [TB^. The upper limit is consistent with theory prediction. 
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Fig. 2. Panel (a): Fit as a function of u^^ijucc- Straiglit lines is for the 
— Xmin + 1- Panel (b): B- (solid curve) /D- (dashed curve) mesons rapidity 
distributions from PYTHIA. 

The upper limit of i? — > e contributions as a function of pr is shown in Fig. [1] 
(b). It is increasing and becomes flat around 7 GeV/c. The pr crossing point, where 
the bottom contribution is equal to charm, of electron spectra from B,D decay is 
very sensitive to the cross section ratio, since at high pT, these electron spectra 
shapes are similar. From the B + D ^ e tit to STAR p+p data, we estimate the 
crossing point p^ > 7 GeV/c. 

Table [5] lists the crossing points of heavy flavor decay electrons in several px 
bins. 



Table 2. Crossing points of heavy flavor decay electrons as a function oi pr- 



PT (GeV/c) 


2 


3 


4 


5 


6 


HPt) 


8 


{B - e)/{D e) < 


0.11 


0.31 


0.53 


0.77 


0.85 


1.2 


1.1 



3. Fit to non-photonic electron V2 

Besides the non-photonic electron spectrum, the non-photonic electron V2 has also 
been measured in 200 GeV Au-|-Au collisions at RHIC [\7\. In this measurement, 
bottom contribution has not been separated, which can be studied by comparing 
simulations and data. Since heavy flavor hadrons pt distributions and V2 are un- 
known, our simulations have to base on the following assumptions: 

— The same relative {B ^ e)/{D ^ e) ratio from p+p to Au-t-Au. 
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Fig. 3. Panel (a): Assumed B- meson V2 (open circles) and D-meson V2 (dashed 
curve) as light meson w^. Panel (b): Electron V2 from B-meson decays (open circles) 
and D-meson decays (open squares). 



— Assume the B-/D- meson V2 as the inputs for the simulation, here we assume 
three aspects: 

• I: B-/D- meson V2 are similar as light meson V2- 

• II; D-meson V2 as light meson V2 but B-meson does not flow. 

• III: B e contribution is neglected and D-meson V2 decreases at pr > 2 
GeV/c. 

Here heavy flavor baryons, Ac, A^ are taken into account as 10% of total heavy 
flavor hadrons [ I17[ ?] . Their V2 are assumed to follow light baryon V2 ■ This baryon 
contribution effect in this simulation is small. 

We use the light meson V2 curve from fitting experimental data [ [19] as the 
input B/D V2 distributions (Assumption I), see Fig. [3] (a). That means in each px 
bin, the B/D A0 distribution is initialized. The electron A(f> distributions in each 
Pt bin will be obtained via B/D decays in PYTHIA model. Then the electron V2, 
shown in Fig.[3](b), will be extracted by fitting the A(f> distributions in each px bin. 

Fig. [3] shows the obvious mass effect: The B/D V2 are assumed as the same, 
but the decay electron V2 can be very different due to decay kinematics. This is not 
surprising, since we know B-meson is much heavier than D-meson and light hadrons. 
The decay electrons can only have a small momentum fraction of B-mesons. The 
momentum and angular correlations between decay electrons and B-mesons are 
weak, especially at low px- Therefore, at low px the decay electron angle will 
almost randomly distribute. So we see the zero or negative V2 for the electron from 
B-meson decays. But from previous study, we know that bottom contribution below 
3 GeV/c is small, thus the mass effect to the total electron V2 is not significant. 

Fig. [4] (a) shows the total electron V2 from PYTHIA simulation compared to 
data. The measured non-photonic electron V2 from PHENIX is shown as the tri- 
angles. The solid curve (Assumption I) is the sum V2 of the two decay electron 
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Fig. 4. Panel (a): The total electron V2 from PYTHIA simulation assuming that 
bottom flows (solid curve) and bottom does not flow (band) compared to data. 
Panel (b): The total electron V2 from PYTHIA simulation fit to data and the 
estimated D- meson V2- 

V2 distributions in Fig. [3] (b) by taking the relative ratio of {B e)/(Z) e) 
into account. It can not describe the data. If we assume B-meson does not flow 
(Assumption II), the total decay electron V2 will become decreasing, shown as the 
band. The band is corresponding to the a^^i/acs = (0.3 — 0.7)% (The upper limit, 
0.49%, is in between). It has better agreement with data, but still higher. The 
decreasing of non-photonic electron V2 could be due to -B ^ e contribution and 
B-meson V2 could be very small. But below 3 GeV/c, B ~* e contribution is not 
significant. That indicates D-meson V2 should be smaller than light meson V2 and 
start decreasing at higher pT (> 2 GeV/c). 

So ignoring _B — > e contribution, we try to speculate the D-meson V2 by fitting 
the data using decay electron V2 (Assumption III). In Fig. |4] (b), the best fit of the 
decay electron V2 is shown as the open circles. The estimated D-meson V2 is shown 
as the dashed curve, which is smaller than light meson V2 above 1 GeV/c and start 
decreasing above 2 GeV/c. 

4. Conclusions 

Charm/bottom and their decayed electron spectra and V2 have been studied us- 
ing PYTHIA simulation. From fitting to the STAR non-photonic electron spectra 
in p+p collisions, we estimate the upper limit of the total cross-section ratio as 
o'fcb/o'cc < (0.49 ± 0.09 ± 0.09)%. And the crossing point of electron spectra from B 
decay and D decay is estimated as > 7 GeV/c. 

The bottom contribution due to mass effect can decrease the non-photonic 
electron V2, but this effect is not significant. The decrease of the non-photonic 
electron V2 is mainly due to the decrease of the parent D-meson V2- The estimated 
D-meson V2 is smaller than light meson V2 above 1 GeV/c and start decreasing above 
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2 GeV/c. Tliis most possible D-meson V2 distribution sliows tliat at < 3 GeV/c, 
wfiere tlie bottom contribution is negligible, D-meson has large V2, indicating that 
charm strongly flows in high dense medium, which could be the evidence of light 
flavor thermalization in QGP created at RHIC energy. 
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